Perfect Absorption in Ultrathin Epsilon-Near-Zero Metamaterials Induced by Fast- Wave 

Non- Radiative Modes 



O 
(N 

o 

Q 



Simin Feng' * and Klaus Halterman' 

^Michelson Lab, Physics Division, Naval Air Warfare Center, China Lake, California 93555 

(Dated: December 6, 2011) 

Above-light-line surface plasmon polaritons can arise at the interface between a metal and e-near-zero meta- 
material. This unique feature induces unusual fast-wave non-radiative modes in a e-near-zero material/metal 
bilayer. Excitation of this peculiar mode leads to wide-angle perfect absorption in low-loss ultrathin metama- 
terials. The ratio of the perfect absorption wavelength to the thickness of the e-near-zero metamaterial can be 
as high as 10*; the electromagnetic energy can be confined in a layer as thin as /!/ 10000. Unlike conventional 
fast-wave leaky modes, these fast-wave non-radiative modes have quasi-static capacitive features that naturally 
match with the space-wave field, and thus are easily accessible from free space. The perfect absorption wave- 
length can be tuned from mid- to far-infrared by tuning the e » wavelength while keeping the thickness of the 
structure unchanged. 

PACS numbers: 42.25.Bs, 78.67.Pt, 42.82.Et 
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A metamaterial is a composite structure with an electro- 
magnetic (EM) response not readily observed in naturally oc- 
curring materials. Many remarkable phenomenon have been 
predicted [1^] by tuning the permittivity e and permeability fi 
in extraordinary ways. If the dielectric response is made van- 
ishingly small, creating an epsilon-near-zero (ENZ) material, 
interesting radiative effects are expected to occur [5-7]. On 
the other hand, by manipulating the EM response to achieve 
a small transmittance (T) and reflectance (R), enhanced ab- 
sorption can ensue. Strong absorption in a thin layer typically 
requires high loss. One of the earliest absorbers, the Salisbury 
screen [8], is based on the phenomenon of destructive wave in- 
terference, and is thus limited to a minimum thickness of one 
quarter wavelength. To overcome the thickness constraint, ab- 
sorbing screens using metamaterials [9] and high impedance 
ground planes [10] have recently been proposed. By explor- 
ing resonant enhancement, thin metamaterial and nanoplas- 
monic absorbers were demonstrated in structures having lo- 
calized resonances [11-23]. In those structures, the geometri- 
cal quality factor (GQF), i.e. the ratio of the perfect absorption 
wavelength to the thickness of the medium, is significantly im- 
proved compared to the standard Salisbury screen. The best 
GQF of those structures is about 40 [14]. 

Based on a fundamentally different mechanism, in this pa- 
per we demonstrate wide-angle perfect absorption in a low- 
loss ultrathin ENZ-metamaterial/metal bilayer as shown in 
Fig. 1. In our structure, the GQF can be as high as lO'*. This 
structure possesses fast-wave non-radiative (FWNR) modes 
due to unconventional above-light-line surface plasmon po- 
laritons (ALL-SPPs) at the ENZ-metal interface, which is eas- 
ily accessible from free space. Fast waves have phase veloci- 
ties exceeding the speed of light in vacuum. Conventional fast 
waves in planar structures are radiative leaky modes that can- 
not be excited by plane wave incidence due to wave-vector 
mismatch at the boundary. For our ENZ-metal structure, 
the exotic FWNR mode is naturally matched to free space 
at the ENZ-air interface, and thus is easily accessible from 
free space. Unlike the Salisbury screen where the thickness 



of the absorbers increases with the absorption wavelength, in 
our structure the wavelength can be tuned from mid- to far- 
infrared (IR) by tuning the e x Q wavelength while maintain- 
ing the thickness. 
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FIG. 1. (Color online) A plane wave is incident (at an angle 9) on an 
anisotropic (ei , = 1 and ei; a 0) ENZ/metal bilayer with thicknesses 
di and A for the ENZ medium and metal (Ag), respectively. Above- 
light-line surface plasmon polaritons (ALL-SPPs) are excited at the 
ENZ-metal interface. 

To investigate this absorption phenomenon in a general 
fashion, an anisotropic ENZ metamaterial (ei t - I, ei~ ~ 0) is 
assumed (see Fig. 1). In the following, the subscripts 1 and 2 
refer, respectively, to the ENZ medium and metal. The areas 
above the ENZ material and below the metal are free space 
and refer to regions and 3, respectively. Since the ei^ has no 
significant impact on the result, eix = I is assumed through- 
out the paper The permittivity of silver in the infrared re- 
gion is obtained by curve fitting experimental data [24] with a 
Drude model, e,„ = 1 - ai-jjap', where the "plasma frequency" 
a>p = 5.38 jum"'. The propagation constant of SPPs at the 
interface of two semi-infinite anisotropic media is given by 
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^ei;yUlv(l -6), 



(1) 



where A:o - w/cand^ = ^ixl^ixK^ixl^u)- The approximation is 
taken when ei - — » 0. In our case, //i ,. - fi2y = 1 andea^ = eiz - 

e,„ lO*" (in the IR region). Thus, kspp/ko ^ yje^^l^ < 1, 

which characterizes above light-line SPP dispersion, and is 
clearly different from conventional SPP dispersion. 
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The absorption can be calculated from Maxwell's equa- 
tions. Assuming a harmonic time dependence exp(-/wr) for 
the EM field, we have 



V X £) = klil, 



E), 



VxH)= klifin ■ H) , 



(2) 



where e„ and Ji„ are, respectively, the permittivity and perme- 
ability tensors for a given (uniform) region (n = 0, 1,2, ■ ■ 
which in the principal coordinates are described by. 



e„ = en.cxx + enyyy + e„-zz , 



(3) 



We consider TM modes, corresponding to non-zero field com- 
ponents Hy, Ey, and E~. The magnetic field Hy satisfies the 
following wave equation: 

1 d^Hy 1 d^Hy - 

which admits solutions of the form iff{z) exp{i/3x). The paral- 
lel wave-vector /? is determined by the incident wave, and is 
conserved across the interface. 
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where a,, is the wave number in the z direction. The func- 
tional form of if/iz) is either a simple exponential exp(/Q'„z) 
for the semi-infinite region or a superposition of cos(q'„z) 
and s,m{anz) terms for the bounded region (along the z direc- 
tion). The other two components E-c and E~ are found using 
Maxwell's equations. By matching boundary conditions at the 
interface, i.e., continuity of //, and E^, the transmittance {T) 
and reflectance {R) can be calculated via the Poynting vector 
S, given by S = c/87r!R(£' x H*). The fraction of energy that 
is absorbed by the system is determined by the absorptance 
(A), where A - I - T - R, consistent with energy conserva- 
tion. Figure 2 shows the absorptance of the ENZ-metal struc- 
ture versus the angle of incidence (AOI) in the mid- to far-IR 
regime. Unless stated otherwise, the results that follow corre- 
spond to ei; = 0.001 + /O.Ol . The AOI at which the perfect ab- 
sorption occurs depends on the wavelength and the thickness 
of the ENZ layer. The full-width-half-maximum (FWHM) an- 
gular width can be as high as ~ 60° (top-left panel), while the 
GQF can reach ~ 10^ (top-right panel). When the thickness 
increases, the angular bandwidth varies differently in the mid- 
and far-IR regions. 

In Fig. 3, we show how the absorptance varies with the 
AOI and incident wavelength (right panels). The left set of 
panels corresponds to the incident angle Op that results in per- 
fect absorption, as a function of wavelength. Perfect absorp- 
tion occurs when the dependent effective impedance iXe) 
of the ENZ-metal structure matches that of free space iZa), 
i.e. Xe - Z,{), in which case the reflection coefficient is zero. 
The effective impedance of the ENZ-metal structure can be 
expressed as 

1 - ri2 exp(/20) 



Z.e - Zl 



1 -I- ri2 exp(/20) 
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(6) 
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FIG. 2. (Color online) Absorptance vs. angle of incidence (see 
Fig. 1) at the wavelength A = 10.94 /jm (left panels) and A = 
201.66;/m (right panels) for two different thicknesses of the ENZ 
medium: di = 0.02 jim (top panels) and di = 0.2 yum (bottom pan- 
els). The Ag-layer thickness (di) is 20 nm. 
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FIG. 3. (Color online) Absorptance (right panels) and direction 
of perfect absorption (left panels) when the ENZ-thickness di = 
0.02 yum (top panels) and rfj = 0.2/jm (bottom panels). Color-bars 
represent the magnitude of the absorptance. The curves in the left 
panels are computed using three different methods (see text). Solid 
(blue): numerically extracted from the corresponding 2D plot in the 
right panels. Dashed (green): obtained from Eq. (6). Circles (red): 
computed from Eq. (7). The Ag thickness is 20 nm. 



where m = (Zi - Z2)/(Zi + Z2), Zj = aj/ej, (j = 0, 1,2), 
and (p - aidi. Perfect absorption occurs at the particular AOI 
that satisfy Eq. (6). It is also informative to analyze the cor- 
responding waveguide modes, which are solutions to the tran- 
scendental equation. 



tanffl = 



. Zi(Z() + Z2) 
Z? + Z0Z2 



(7) 



which implicitly depends on the parallel complex propaga- 
tion constant /3 (Eq. (5)). The mode solutions of Eq. (7) have 
four branches due to the + signs inherent to the square root 
of ffo and 02 in Eq. (5). The branch leading to perfect ab- 
sorption corresponds to both ao < and 02 < 0. In the long 
wavelength regime and for ultrathin slabs, this branch pro- 
vides fast- wave (fi/ko < 1), non-radiative (J3 real) modes that 
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are guided along the ENZ layer. In the left panels of Fig. 3 
we compute the direction of perfect absorption as a function 
of wavelength in three complimentary ways: The solid (blue) 
curves are numerically extracted from the absorptance corre- 
sponding to the right panels, the dashed (green) curves are ob- 
tained from Eq. (6), where the structure is impedance matched 
to free space, and the circles (red) are calculated using Eq. (7) 
and a root searching algorithm. The overlap of each set of data 
further demonstrates that the underlying mechanism of perfect 
absorption is the excitation of FWNR modes. We see also that 
perfect absorption occurs only at oblique incidence where the 
electric field has a nonzero component in the normal direction. 

The energy density, U, in lossy and dispersive anisotropic 
media is given by. 



(8) 



We show in Fig. 4 the spatial distribution of the energy den 
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FIG. 4. (Color online) Spatial distribution of the energy density com- 
puted from Eq. (8) (normalized so that Hy is unity at the ENZ-air 
interface) for X = lOyum and AOI = 16.1° (left) and A = 155 /jm and 
AOI = 56.3° (right). The ENZ-layer with d, = 0.2 jim is located at 
z e [-0.1 0.1]. The regions z e [0.1 0.15] and z e [-0.2 - 0.1] 
are, respectively, silver and air. Color-bars represent the magnitude 
(xO.Ol) of the energy density, indicating the energy is more confined 
to the ENZ medium at far-IR wavelengths. 

sity in the ENZ-metal structure under conditions of perfect 
absorption at mid- and far-IR wavelengths. In the mid-IR (left 
panel), the field inside the ENZ-layer is concentrated at the 
two interfaces, and a very thin layer of ALL-SPPs can be seen 
at the ENZ-metal interface. The strong field at the ENZ-air 
interface is distinguishable from conventional surface waves, 
since in this case the exponential decay of the field occurs 
strictly inside the ENZ medium. Outside of the ENZ medium, 
the field neither decays like a surface wave nor does it in- 
crease like a leaky wave. The observed uniformity indicates 
that the field outside the structure is plane wave matching with 
the free-space wave, and thus can be easily excited by plane 
wave incidence - a striking difference from surface waves and 
leaky waves. In the far-IR (right panel), the field is even more 
uniform and confined to the ENZ-medium and also matched 
to the space-wave field at the ENZ-air boundary. Around the 
point where the real part of e approaches zero, even with a 
small loss 3(ei;) - 0.01, the presence of the FWNR mode 
reinforces the increasing absorption due to the strong field in- 
side the ENZ medium (by virtue of the electric displacement 
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FIG. 5. (Color online) Direction (top left) and FWHM bandwidth 
(bottom left) of perfect absorption vs. ENZ thickness. The curves are 
numerically extracted from the absorptance when A = \Ojim (dots 
(blue)), A = 100 yum (solid (green)), and A = 300 /jm (dashed (red)). 
The circles in the left- top panel are obtained from Eq. (7). The nor- 
malized propagation constant p/ko corresponding to these modes is 
shown in the right panels. Top-right: real part of /3/ko. Bottom-right: 
imaginary part of p/kg. The Ag thickness (A) is 10 nm. 



continuity) leading to perfect absorption in the ultrathin ENZ 
medium. 

The minimum thickness of ENZ materials is generally lim- 
ited by fabrication methods except in some cases of naturally 
occuiTing resonances. Figure 5 shows the direction, 9p (top 
left), and FWHM angular bandwidth (bottom left) of perfect 
absorption versus ENZ thickness for three wavelengths, as 
well as the coiTesponding normalized complex propagation 
constants (right panels). The circles in the top left panel repre- 
sent solutions to Eq. (7), which are consistent with the absorp- 
tion calculated using the Poynting vector These curves also 
correlate with the solutions from the matched impedance ex- 
pression (Eq. (6)) (not shown). When the thickness increases, 
the direction of perfect absorption is shifted towards the sur- 
face normal for both mid- and far-IR regions. There is an 
optimal thickness that yields the largest bandwidth. The real 
part of /3/ko is above the vacuum light-line, indicating a fast- 
wave character for all three wavelengths shown. The imagi- 
nary component of /3 is zero for the far-IR wavelengths. For 
^ = 10pm, 3(j6) = when di < 0.2pm. When di > 0.2pm, 
f3 starts to pick up a small imaginary part which increases with 
thickness: The fast-wave non-radiative mode transforms into 
a fast-wave leaky mode [7] as the thickness increases. 

To further understand the FWNR modes, we integrated the 
transport power S x, dissipation rate D, and %{E,) over the 
ENZ-thickness (normaUzed by the input £;). The dissipation 
rate is given by. 
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Stt L 



■fiH]. 



(9) 



The average energy density in the ENZ medium is also calcu- 
lated (Eq. (8)), as well as the effective surface charge density 
Q - €i,E,/4n (normalized by the input E~) and capacitance 
C - Q/V per unit area, which is much larger than the capaci- 
tance (Co - e\z/di) of a parallel-plate capacitor (of unit-area). 
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FIG. 6. (Color online) Average energy density (left top), transport 
energy (5 , /c) per unit area (left bottom), voltage (V) across the ENZ- 
layer (middle top), the ratio of power to dissipation rate (middle bot- 
tom), the effective surface charge density Q (right top) and the nor- 
malized capacitance C/Co (right bottom) vs. d[ for A = lOyum (dots 
(blue)), A = 100 yum (solid (green)), and A = 300 /jm (dashed (red)). 




FIG. 7. (Color online) Color map of the absorptance versus AOI 
and loss (3(ei;)) for A = lOyum (top panels) and A = 200 /jm (bottom 
panels). The left and right set of panels correspond to ^R-(ei;) = 0.001 
and ^!l(ei;) = 0.01 respectively. Here, di = 0.15 yum and d2 = 10 nm. 
The color bars represent the magnitude of the absorptance. 

due in part to the field enhancement inside the ENZ mate- 
rial. These results are presented in Fig. 6. The confinement 
is stronger for far-IR and decreases with the thickness (left- 
top). The linear responses of the voltage (middle-top) and 
constant behavior of the charge (right-top) as a function of 
thickness are hallmarks of capacitors. The deviation from the 
linear behavior for A - 10 /im can be understood as a balance 
between radiative loss, which acquires more "charge" stored 
in the "capacitor" (right- and middle-top) and energy transport 
(left-bottom) (see right-bottom panel of Fig. 5 where 3(J3) + 
when d\ > 0.2 yum for A - lOjum). The radiative loss reduces 
the enhancement of the capacitance at A = lOjum as shown 
in the right-bottom panel. In essence, the higher the confine- 
ment, the larger the capacitance. The effective capacitance 



is found to be equal to cQ^ l{2S x) or 25.j/(cy^), which im- 
plies the transport energy is driven by quasi-static "surface 
charges". Thus, FWNR modes can be considered quasi-static, 
capacitively driven fast waves. The middle-bottom panel of 
Fig. 6 shows the ratio of transport power to dissipation power, 
suggesting longer wavelengths have better power handling ca- 
pability if the structure is used as an ultrathin channel to trans- 
port light. The influence of loss on the absorption is illustrated 
in Fig. 7 for different ei; (real part) in the mid- and far-IR. In 
the far-IR, the smaller V>(€\.) has larger angular bandwidth, 
which is opposite to what occurs in the mid-IR. 

In summary, wide-angle perfect absorbers have been 
demonstrated in ultrathin ENZ-metal structures, and which 
have extraordinary high geometric quality factors. Such struc- 
tures may have applications involving ultra light-weight in- 
frared absorbers, cloaking, EM shielding, photovoltaic sys- 
tems, and highly efficient infrared sensors and detectors. The 
authors gratefully acknowledge the sponsorship of ONR N- 
STAR and NAVAIR's ILIR programs. 
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